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1. Introduction 



Given a quantum state \ y) and a corresponding density matrix p = \ l^t us introduce 

two observers A and B = A, such that observer B sees only degrees of freedom complimentary to 
degrees of freedom seen by observer A. Then one can form a reduced density matrix for observer 
A: Pa = Tr^p, and define the entanglement entropy as the von Neumann entropy of reduced states 



Sa 



-TvAPA^ogpA = -J^^dogXi, 



(1.1) 



for a 



where A,- are eigenvalues of the reduced density matrix. The entanglement entropy Sa 
product state, and is at maximum for a maximally entangled state: < 5a < 1 [ebit]. 

As a simple example consider a bipartite system = cos 0| Ia) ® I Ib) + sin 0| Ia) | Tb). 
The reduced density matrix is Pa = cos^ 0| |)(| | + sin^ 0| |)(| |, while the entanglement entropy 
is 

5a = — 2cos 01ogcos — 2sin 01ogsin0, (1.2) 

and takes its maximum value of log 2 when cos^ 6 = j. 

One of the more interesting application of the entangle- 
ment entropy is its use as a probe of phase, which is espe- 
cially important in situations where an order parameter is not 
known. It turns out that in confining gauge theories the en- 
tanglement entropy may show an interesting non-trivial be- 
havior. First it was studied in gravity duals of confining large 
gauge theories [|l|], where the following geometry was con- 
sidered (see Fig. |]) 

A 

B = M^-i X (M-I,), (1.3) 

here 1/ is a line segment of length I. It was observed that 
at length separation / = 1^ the entanglement entropy exhibits a non-analytical change in behavior, 
reminiscent of a phase transition. 

A natural question to ask if this behavior is relevant for 
small N gauge theories. This question was first addressed 
in Migdal-Kadanoff formalism |^. Subsequently other mea- 
sures of entanglement (purity) for SU (2) gauge theory were 
studied in Monte Carlo simulations [^. The results of the 
Migdal-Kadanoof study are reviewed in this proceedings. 

The expression for the von Neumann entanglement en- 
tropy is obtained using the replica trick [Q, ^]: « replicas of 
the system (with region B integrated out) are glued along the 
time boundaries of region A. We demonstrate the procedure 
for a 2d system in Fig. ^ Taking the trace of the combined 
system we obtain its normalized partition function 

Z„(A) 



Figure 1: Partitioning by two imag- 
inary surfaces into regions A and B. 



R 



T 



1/T 



Trp^ 



Z" 



(1.4) 



Figure 2: Z„ for 1 + 1 dimensional 
gauge theory. 
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where Z = Z\ is the partition function of the original system. Then after performing analytical 
continuation to real n one obtains the entanglement entropy 



d „ d ZJA) 
Sa = - lim ^Trp; = - lim ^ 



(1.5) 



We apply the described technique to SU (N) gauge theory inD = d + l dimensions 

Z= [YldU,Y[e-'', (1.6) 

I p 

where Sp = S{Up) = —^/{2N)TrUp + h.c. is a standard Wilson lattice action with the plaquette 
variable Up = Wi^dp^i ^^id ^ = '^N /g^ is the inverse lattice coupling. The gauge invariant action 
is a class function and therefore allows for character expansion 



r \ r^O I 



(1.7) 



where c, = F./Fq < 1 and F, = j dUe-^^^^}-^x;{U). 

2. SU{N) gauge theory in 1 + 1 dimensions: Exact solution 

We start with exactly solvable 2-dimensional SU {N) gauge theory (for an overview and large 
N treatment of zero temperature U (N) gauge theory see Ref. [^). At finite temperature the gauge 
theory leaves on an M x Si manifold compactified in time direction with period l/T. The corre- 
sponding lattice theory is formulated on an Nr x Nt time direction periodic lattice, with space-time 
cut-off a and aNf = l/T and aNr = R. 



After the partition function is character expanded ( |1.7| ) one can consider various contributions 
to it from minimal surface elements bounded by a single loop dA 



fi{ay,dA) = l + Y^aidiXi{dA), 



(2.1) 



The contribution to the partition function of two such surface elements A and B with a common 
boundary AHB, which is integrated out, defines a new elementary surface 



f{{c};d{AUB)) = j d{Ar\B)f{{a]-dA)fm\dB) = 1 + £c,-J,x(5(AUB)), (2.2) 

where c, = aibi. Thus the junction of the surfaces in the space of character coefficients is repre- 
sented by an ordinary product. 

Now for any 2-dimensional surface one can integrate all the internal links (in this way joining 
elementary surfaces). The resulting expression for the partition function is 



Y[dUiY,FXXriUdA), 

led A r 



(2.3) 



3 



Entanglement entropy in SU(N) gauge theory 



Alexander Velytsky 



where A = NrNt is the area of the total surface in plaquette units and dA is the contour enclosing 
the surface. A similar result follows for the partition function of the glued system Z„, but with the 
corresponding surface area A„ = nA = nNrNt and perimeter 5A„. 

In order to evaluate the perimeter integral one needs to set the spatial boundary conditions 
(BC). The free BC produces a trivial result Z = F^ and Sa = 0. Therefore we set the periodic BC. 
The perimeter integrals for Z and Z„ result in 



Z : jdV j dUXriUVU'^V'^) = [ dV^Xr{V)Xr{V'^) = ^, 

J J J dr dr 

Z,:J^rj ^rj ^ X>-{Ul)...Xr{Un)Xr{Uj)...Xrm) _ 1 



dU\...dUn 



dr 



n—l 



(2.4) 
(2.5) 



and one can obtain the entanglement entropy 



d_Zn 

dnZ" 



n=\ 



\og{l + Y,Cr) ^.y^,A ■ 



(2.6) 



It is interesting to note that 5^(0 is /-independent |/ / 0, and its behavior is similar to end-point 
phase transition Sa (/ = 0) = 0. 

In the expression for the entanglement 



entropy ( |2.6| ) one can truncate the series if 
A > > 1 or, equivalently, if j8 is small. Let us 
be more specific and consider SU (2) gauge 
theory. Then Cr = /2r+i(2j8)//i(2j8), where 
In{x) is the modified Bessel function. In 
strong coupling limit 



1 



j2r+l 



/2r+l(2iS) 

where < 2j8 << V2r + 2, 



r(2r + 2)' 



(2.7) 



and Cr = j8^'^/(2r+ 1)!. One can then sum 
terms of expression ( [2.6[ ) to a given precis- 
sion to get 




Figure 3: SU{2) entanglement entropy at various lat- 
tice j3 and surface area A. 



Sa = £(l+log4-log£), e={hr 



(2.8) 



We present the results for the entanglement entropy of 1 + 1 dimensional SU (2) gauge theory 
as a function of the lattice coupling for various surface areas A = 1,2 and 10 in Fig. ^. For A = 1 
we also plot the strong coupling perturbation result. 
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2.1 Higher dimensions: Migdal-Kadanoff treatment 



The higher dimensional SU (N) gauge theories cannot be 
solved exactly, however the Migdal-Kadanoff (MK) approximate 
approach is known to produce robust results in studies of phase 
structure of the models. The standard MK decimation procedure 
(A -transformation) in D = d + 1 dimensions is defined by the 
following recursive steps (see Fig. ^) 



. r 

Fr = jdUe-^"'"(^)jXm- 



(2.9) 




Figure 4: Illustration of the 
Migdal-Kadanoff procedure. 



Here A is the scaling factor of the renormalization group (RG) 
transformation, = A^^^ is the factor by which we strengthen 
the interaction on the resulting coarser lattice and A = A^ is the 
surface of the new elementary plaquette in units of plaquettes of 

the underlying fine lattice. One can choose between two orderings of bond moving and strength- 
ening of the couplings steps: b = corresponds to Migdal, while b = \to Kadanoff prescription. 

As an example consider a gauge theory formulated in a 2 + 1 dimensional box. One can then 
decimate out all bulk degrees of freedom and be left only with the perimeter integral. Setting free 
boundary conditions one then gets 



Z = f dUdVf{{c^,iy,U^VUV^)f{{ct,i];V) = 1 + Y,c^,i + £ c^AdjCt,p\j, (2. 



10) 



where 



f{{c,}-dA) = l + Y.diCr,iXi{dA, 



(2.11) 



and 



4 = JdVXk{V'^)XiiV)XjiV) 




are coefficients of the Clebsch-Gordan series 
of irreducible representations (see for details). 

Now we can apply MK procedure in order to com- 
pute the partition functions for ordinary Z and glued Z„ 
systems, see Fig. In order to cancel out the contribu- 
tion from the bulk we cary out decimations for Z„ and 
Z in exactly the same manner. First we start with the 
standard MK decimation procedure (A -transformation) 



(2.12) 

Y.k^ii^^^'^ for the Kronecker product 




(2.S). The decimation should be altered when the lat- 
tice spacing becomes equal to / (the smallest scale in 
the problem). This affects only /-plaquettes, which have 





















1 









1/T 



R 



Figure 5: Z„ for 2 + 1 dimensional theory. 
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one link spanning the entangled region. In this case one 
still can move plaquettes in D — 2 direction but the tiling 
is done with A plaquettes 

r 1 

l^F.^drXriU) , Fr = JdUe-^'''':'<^''^jX:{U). (2.13) 



We refer to such a procedure as p -transformation. All the other plaquettes are unaffected by this 
change and are decimated according to the standard (A -transformation) procedure. 

If we choose / to extend in y direction, we can write the final formula for the entanglement 
entropy [§] 

SA = -P,,o + logZ-^ (2.14) 



where the dot stands for X = 4-X 

dn 



and 

n=l 



fn = I«,cl,C3,,,f (log Ifk^^ [ 1 + £ djCrjDij] + £ djtrjDij] . (2.17) 

The final expression appears to be very complicated. However, one can analyze its behavior as 
a function of I. First we simplify the expression by choosing the geometry with specific symmetry, 
so that c*; = ; = Note that the dependance on I is encoded in the values of c-. Note that / 
regulates when A -transformation is switched to p -transformation, i.e. it sets the initial value for 
c^(mo) under p -transformations. Next we analyze the RG flow of SU (2) gauge theory for c- (m) as 
a function of number of iterations m under Migdal recursion depending on the starting point. We 
simplify the numerical simulation by considering a starting action in the wilsonian form on Ni^t = 1 
lattice. Note however that even after the first iteration step the action is of a single plaquette form, 
but generally has all irreducible representations. 

We show the projection of the RG flow on the fundamental-adjoint plain in Fig. ^. We observe 
that depending on the starting coupling j3, which sets the /-scale, the flow is in different directions. 
This is a clear indication of a phase transition like behavior. For the scale factor we choose A = 1 . 1 , 
which is known to map correctly the mixed action phase diagram [Q]. The transition occurs at 
j8c G 0.(62,0.65), which corresponds to length scale /* 

Z,7Zc G (1.56,1.66), (2.18) 

where Zc = 1 /Tc is the QCD scale. 



3. Summary 

We studied the entanglement entropy in <i + 1 SU (N) gauge theory. The d = I theory was 
solved exactly. Setting periodic BC we obtained a non-zero universal value for the entanglement 
entropy, which is independent of the size / (end-point transition). 
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Figure 6: Migdal decimation flow for 3 + 1 dimensional SU{2) gauge theory. Projection to cj^j ^"d c\; 
()3 , A ) are indicated. 

Using MK decimation we approximately computed the ratio of partition functions and en- 
tanglement entropy for d > 2. For 3 + 1 SU (2) gauge theory we demonstrated that there is a 
non-analytical change in the RG flow of character coefficients c which define Sa- This allowed us 
to find the length scale where the transition occurs l*/lc G (1.56, 1.66), which is comparable to the 
value l*/lc = 2 derived for infinite A'^ theory 

It is worth to point out that other measures of entanglement, such as the purity jj. = Trp^, 
Tsalhs entropy Sq = l/(<7— 1)(1 — Trp^) and Renyi entropy Hi = I / {I — q) log Trp^ (specifically 
atq = 2) are all dependent on c- and therefore will show non-analyticity in the RG flow. 
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